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Guiding of Whistlers in a Homogeneous Medium 

R. L. Smith 

(May 4, 1960) 

The velocity of energy flow of whistlers in a homogeneous medium is computed as a 
function of wave-normal angles. The maximum allowable cone of ray angles approaches 
19° 29' at very low frequencies, decreases with frequency to a minimum of 11° at a wave 
frequency of one-fifth the gyrofrequency, then increases to 90° at the gyrofrequency. The 
velocity of energy flow departs markedly from the longitudinal value except at very low 
frequencies or very small wave-normal angles. 



1. Introduction 

Part of the energy from a lightning stroke can 
penetrate the lower ionosphere and be guided 
approximately along the earth's magnetic field in 
the outer ionosphere. The energy is dispersed as 
it travels in this region, causing the original signal, 
which can be considered an impulse, to be stretched 
or dispersed into a gliding tone, typically lasting 
about a second. This gliding tone is called a 
whistler. The dispersion is a measure of the 
electron content along the path. 

Storey [1] 2 has analyzed the properties of whistler 
propagation for the case of wave frequency very 
low compared to the local electron plasma frequency 
and gyrofrequency. Briefly, his results show that 

(a) the direction of energy flow lies within a limiting 
cone of 19°29' around the magnetic field direction; 

(b) the longitudinal expression for group velocity 
can be used to describe 1 the velocity of the energy 
flow with an error less than 8 percent for wave 
normal angles up to 70°; (c) if the initial wave 
normal angles for the different frequency components 
are identical, the path of energy flow is independent 
of frequency; and (d) the time delay T from the 
causative lightning stroke is given by T=D /""**, 
where / is the frequency and D is a constant for a 
particular whistler. This relation is sometimes 
called "Eckersley's law". Storey found good agree- 
ment between this relation and data he obtained at 
Cambridge, England. 

From recordings made at high latitudes, a new 
phenomenon, called the nose whistler, was discovered 
by Helliwell et al. [2] in 1955. These whistlers 
exhibited simultaneous rising and falling tones 
starting at the frequency of minimum time delay 
called the "nose". The nose whistler occurs when 
the wave frequency of the whistler becomes com- 
parable to the local gyrofrequency at the top of a 
field line. Even at the lowest frequencies of most 
nose whistlers the departures from Eckersley's law 



i Contribution from Radioscience Laboratory, Stanford University, Stanford, 
Calif. The main results of this paper were presented at the Symposium on VLF 
Radio Waves held in Boulder, Colo., January 1957. 

2 Figures in brackets indicate the literature references at the end of this paper. 



are significant. Other whistlers, though they may 
not exhibit a nose, often show similar departures 
from Eckerslcy's law. These deviations arc im- 
portant in determining the location of the path of 
propagation [3]. Without the high-frequency infor- 
mation, it is usually necessary to assume that 
whistler energy follows a field line terminating at 
the receiver. The diameter of the effective area of 
whistlers has been shown by Storey [1] and Crary 
et al. [4] to be approximately 1,000 km. Consider 
two possible whistler paths terminating at 55° and 
60° geomagnetic latitude, well within the effective 
area of a whistler. The calculated average electron 
density for a whistler of given dispersion will differ 
by 2:1 for these two paths. Clearly the path needs 
to be defined more precisely. Before we use the 
high-frequency information however, we must re- 
examine the equations used to determine the 
dispersion and path. 

The present paper is an extension of Storey's work, 
removing the restriction that the wave frequency be 
small compared to the gyrofrequency. As frequency 
increases, the analysis shows that the limiting cone 
of rays first decreases to 11° at about one-fifth the 
gyrofrequency, then increases to 90° at the gyro- 
frequency. The error in the quasi-longitudinal ap- 
proximation for group velocity increases with 
frequency and wave normal angle. 

2. Refractive Index for Whistlers 

Experimental evidence shows that the attenuation 
of whistlers in the outer ionosphere is very low. 
There are estimates of attenuation of the order of 10 
db and lower over paths which are roughly 10,000 
km long. The attenuation will therefore be ne- 
glected. The square of the wave refractive index for 
the ordinary mode, as defined by Booker, is then 
given from the Appleton-Hartree equation as: 
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6= angle between the wave normal and magnetic 
field, 



/=wave frequency, 

plasma frequency, 
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iV= density of electrons, 
e= charge of the electron, 
m=mass of the electron, 
H= magnetic field strength, 
€o= permittivity of free space. 

If we assume that tan sin 0<<2 (X—1)/Y H , 
i.e., if we assume quasi-longitudinal propagation, 
then eq (1) is simplified to: 
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For most cases of interest in whistler propagation 
it can be assumed that the square of the refractive 
index is large compared to unity. Equation (2) can 
then be further simplified to : 
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where \=llY H =f/f H = normalized wave frequency. 
This relation defines the wave surface. 

Equation (4) shows that for a given value of A, the 
angle between the wave normal and the field direc- 
tion must be less than that given by cos 6 max =^ for 
a propagating wave. The condition of validity of 
the quasi-longitudinal assumption at the maximum 
angle can easily be shown to be : 
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Assuming the wave refractive index is given by eq 



(3), then the group refractive index, from which we 
determine the group velocity, is: 
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The error in the group velocity as derived from eq 
(3) instead of eq (2) is negligible at \=0, 0.75, and 
1.0. The greatest error occurs at X= 0.283 and is 
approximately 
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3. Behavior of the Ray 

The direction of the ray is the direction of construc- 
tive interference of phase fronts for infinitesimal 
changes in the wave normal direction. If the direc- 
tion of the ray with respect to the wave normal is a, 
and a is positive when the wave normal lies between 
the ray and the field direction, then a is given by 
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This relation has been shown by Bremmer [5], Storey 
[1], and others. 

The total angle between the ray direction and the 
field is 0+a. A sketch of (0+«) as a function of 
and parametric in X is shown in figure 1. 

The diagonal line represents the ray direction when 
the wave normal is at its maximum value. The 
greatest negative value of (6+ a) is given by 
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Figure 1. The ray direction as a function of the wave direction. 
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The greatest positive value of (0+a) is determined by- 
setting 



50 -° 

50 -1 ^d0 L 2(cos0-X)J 

^ 3 cos 2 0— 6X cos 0+4X 2 - 1 = 
2 (cos 0-X) 2 + sin 2 



Solving'eq (7) for cos 0, we find 
cos 0=X 
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The maximum positive value is then easily shown 
to be 

tan (fl+.W ^ 1 -^ 2 -^ 372 . (U) 

\ i /max 2 3/2 (l X 2 ) 3/4 

The maximum ray direction for any given normalized 
frequency is given by the larger of the values deter- 
mined from eqs (8) and (11). Figure 2 shows the 
results graphically. 
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Figure 2. Maximum possible ray direction as a function of X. 



The minimum ray cone of 11° occurs at X= 0.189. 
The value at X=0 is found to be 19°29 / in accordance 
with Storey. The half beamwidth lies between these 
two values for all X between and 0.33. 

The group ray refractive index is given by M' = n' 
cos a. The group ray velocity is the velocity of 
propagation of a point on a wave packet limited in 
both length and width where the phase is stationary 
with respect to independent variations of both fre- 
quency and wave normal direction. From eq (6), 



we can obtain 
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and from eq (7) we can obtain 
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Combining (12) and (13), we obtain 
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and ju/^group ray refractive index of a longitudinal 
wave. The factor $> is plotted in figure 3 as a func- 
tion of 0, with X as a parameter. The curves show 
to what extent the group ray refractive index is inde- 
pendent of wave normal direction. 



2.4 




1 1 


1 


1 














- 


2.2 






















- 


2.0 






















- 


1.8 
1.6 




X- 9 


8 


? 


6 


5 


*/ 


37 


■«?/ 




- 




















/ ' 




1.4 


- 




















1.2 




















- 










V 


1.0 






8 


- 
















\° 


.6 


- 
















\ - 


4 


- 
















\- 


2 


- 


1 I 




I 










1 



30° 40° 50° 60° 

WAVE NORMAL ANGLE 9 



70° 80° 90° 



Figure 3. 



Correction factor $ as a function of wave normal 
direction, with X as a parameter. 
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4. Conclusion 

The spread of allowable ray angles approaches 
19°29' as the normalized frequency X=///h- tends 
toward zero. The limiting cone reaches a minimum 
of 11° at about one-fifth of the gyrofrequency, then 
rapidly increases with increasing frequency. The 
energy near and above the nose frequency of whistlers 
(A — 0.25) would then be expected to diverge consid- 
erably in the outer ionosphere unless an additional 
confining mechanism is postulated. 

Furthermore, the group ray refractive index may 
depart markedly from the longitudinal value. Even 
for wave-normal angles less than one-half the limiting 
angle the correction factor may exceed unity by more 
than 10 percent. Evidence for the spread in group 
ray refractive index should be found in the time delay 
spreading of received whistlers. Since this spread is 
observed to be less than 1 percent in most nose 
whistler traces, a mechanism for confining the wave 
normal angles to small values is indicated. Such a 
mechanism might be enhanced columns of ionization 
alined along the magnetic field as discussed by Helli- 
well [3] and Smith et al. [6]. 
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